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On the Mathematical Theory of Two Star-drifts , and on the 
Systematic Motions of Zodiacal Stars. 

By A. S. Eddington, B.A., M.Sc. 

In a former paper,* the distribution of the proper motions of 
the stars was discussed from the point of view of the hypothesis 
of two drifts of stars. The method of analysis consisted in finding 
by trial and error the constants of two drifts which would lead to 
a distribution of proper motions (as regards direction) agreeing as 
nearly as possible with the observed distribution. There are 
evident disadvantages in a method of this kind; the difficulty of 
simultaneously adjusting five disposable constants, so as to obtain 
the best agreement with observation, is considerable; there is 
room for bias in deciding what is the best agreement with observa¬ 
tion ; and the method is inapplicable to the discussion of a series 
of proper motions of stars sparsely distributed over a large region 
of the sky, such as the Bradley proper motions. 

It is the object of the first part of this paper {Mathematical 
Theory ) to obtain from the observational data direct equations to 
determine the constants of the drifts. In the second part the 
theory is applied to the analysis of. the proper motions of the 
Catalogue of Zodiacal Stars ( Astronomical Papers of the American 
Ephemeris, vol. viii., part iii.); this is an example of a series of 
proper motions whioh, while sufficiently numerous to afford a good 
determination of the velocities of the drifts, are much too sparsely 
distributed to be analysed by the method of trial and error. 


Mathematical Theory. 


I have defined a “ drift ” to be a system of stars in which the 
motions of the individual stars relative to the mean of the system 
are haphazard; the whole system has in general a drift-motion 
relative to the Sun. For the sake of precision, the distribution of 
the haphazard or “ peculiar motions ’’ in the ideal drift is taken to be 
according to Maxwell’s law. 

In “Systematic Motions,” p. 36, the law of distribution of the 
proper motions in the different directions was discussed for a 
system of stars forming a drift thus defined; the analysis may be 
briefly recapitulated here. 

A small region of the sky is considered, and motions in the line 
of sight are neglected, so that the problem is two-dimensional. 

Let the number of stars having component peculiar (linear) 
velocities between the limits ( u, v) and (u + du, v + dv) be 


nff e _ hV+v , )dudv 



7T 


in accordance with Maxwell’s law. 

* “The Systematic Motionsof the Stars > 9 Monthly Notices, vol. ixvii. p. 34.. 
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n 1 is the total number of stars of the drift. 

h a constant inversely proportional to the mean peculiar speed 
of the stars. 

Let v 1 be the velocity of the drift. 

0 1 the inclination of the direction of this drift-velocity 

to the axis of x. 

r, 0 the amount and direction of the resultant linear 
velocity of a star ( i.e . drift velocity and peculiar 
velocity compounded). 

p dO the number of stars having proper motions in 
directions inclined to Ox between 0 and 0 + dO. 

Then u 2 + v 2 = r 2 + vf - 2 v x r cos (0 - Of 

dudv = rdrdO 

and hence pdO = r ^-do( e -A 2 { r 2 +*’i 2 _ 2 , ’i rC 0 B (* -e i )}rdr . ( 2 ) 

TTJo 

By means of the substitution 

x — hr — hv l cos (0 — Of 
this expression can easily be reduced to 

p = i + re^J^e-^dxj . . (3) 

where t = liv l cos (0 - Of 

or setting f(f) = % + re r2 J e~^dx 

p = !h,e- h ^f(hv l co8(0-0f) • • . ( 4 ) 


The function / was tabulated in “ Systematic Motions,” p. 36 ; a 
2 

factor —>- was, however, inserted in defining it there, which it has 
\Jir 

been more convenient to drop now. 

Consider now a series of proper motions of stars belonging to 
two drifts; distinguishing the corresponding constants by suffixes 
1 and 2 , we shall have 


p = ^ e h2v 'V(hvi c os 0 - 0 f + ^e h ^f(hv 2 cos 0 - Of . ( 5 ) 

If the observed proper motions are considered to correspond to 
two drifts, a series of equations of condition will be obtained by 
giving various values to 0 in ( 5 ), the corresponding values of p 
being derived from observation. From these the six unknowns, 
n v n 2 , v v v 2 , 0 V 0 2 are to be derived; these, however, may be at 
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once reduced to five unknowns, for n l + n^ is equal to the total 
number of stars. 

It is of course impossible, with non-linear equations of con¬ 
dition such as these, rigorously to determine the values of the 
constants which shall make the sum of the squares of the residuals 
a minimum. In order to obtain what may serve, as it were, for 
normal equations, I have expanded each side of equation (5) in a 
Fourier series, and have equated the first few coefficients on each 
side until sufficient equations were thus formed. It is clear that, 
at least if the observed distribution really does correspond to two 
drifts, a good agreement between the observed and theoretical 
distributions will have been obtained when the Fourier series 
representing'them agree as far as the coefficients of cos 3# and 
sin 3#, and it is not likely that the agreement could be much further 
improved. Discordances represented by differences in the co¬ 
efficients of higher harmonics are likely to be of an accidental and 
unimportant character. 

In practice, the process of. equating Fourier coefficients is 
equivalent to equating the sums 2 cos 0, S sin 0, 2! cos 20, 2 sin 20, 
etc., for the observed proper motions, to the corresponding integrals 
for the theoretical distributions. These integrals are evaluated in 
the note at the end of this paper. We find there that for a single 
drift (n v v 1 0 X ), 


J™p cos (0 - 0 1 )d0J | J 0 (i$h\ 2 ) - V) } 

= »i 1 0 ( 7 ^) | say * 

C 27T / T — (>— 

J p cos 2(0 — O^jdO — - ^2 ) ==n i\ D(H) 1 8 a y. 

j 2K p cos 3(0 - OJcW = n x 1 J 0 (i^V) ~ tJ x (c^h\ 2 )(i 

= w 1 |E( 7 tw ] )| say 


( 6 ) 



and evidently from the symmetry of the drift 

*2ir f 27T f 27r 

p sin(00 x )cZ0 = I psin 2(6 - 0j)c£0 = I psin 3(0 - 0 1 )d0 = o . . (7) 

o J O J o 

A table of the functions here introduced is given later (see 
Table I.). 

It is here convenient to introduce complex quantities into the 
analysis, which will be denoted always by capital letters. 

Let Yj = v 1 e t0 i 

C 1 =C(7 i V 1 ) = |C(H)|e^ ) , 

D 1 =D(fcY 1 ) = |D(Av 1 )|e a ‘«i V . . . (8) 

E x = E(/iV 1 ) = | E(hv x ) | e3 t(? i j 

The moduli have already been defined by equations (6). 
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Still considering one drift only 

/*27T I" 27T 

j pe‘°d 9 — e t0 i I pe^ ~ e ^d0 

Jo Jo 

f 27T f 2TT 

= e tfl i J p cos (0 - Q-^)dO + ie^i j p sin (6 — 0 x )a!0 
= e t0 m x | C(A« 1 ) |. 

f 27 r 

Thus I pe l9 d 0 = » 1 C 1 

J o 

f 27 r 

Similarly. I pe 2t W = w x D x . . . . (9) 

^ O 
/* 27T 

I pe3 l9 d9 = re x E x 

J o 

Keturn now to the consideration of two drifts. Knowing for 
every star the apparent direction 6 (measured from an arbitrary 
initial line) of its observed motion, we can calculate 2 cos 0, 
5 sin 6, 2 cos 26, etc. for all the stars. 

Let L = — 2 e 10 = — (2 cos 0 +sin 

n n \ J 

M = — 2 e at0 = - (2 cos 26 +sin 2<A 
n n\ J 

N = - = — (iZ cos 3$ +sin 3#^ 

where n is the total number of stars discussed. 


And let li=A(i +a ) 
n 

so that -2 = i( 1 — a) since n = n. + n „. 

n 

The equation of the observed and theoretical values of 2 e t0 , 
Se 210 , and 2e3 t0 , in the case of two drifts, therefore leads to 

2 L = (I + a) C x + (I — a) Cg ) 

2M = (i+ a)D 1 + (i-a)D 2 > . . . (10) 

2N = (1 + a)E x + (i — a)E 2 ) 

To solve these equations write 

C= £EM 



Since the arguments of C 15 D x , and E x are 6 V 2 0 V and 3# x , £ and 
y are real, and the argument of P x is 9 V Further, it is found that 


ft 
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for the range of values of hV (up to about [ AY | = i*8) occurring in 
the actual drifts, y may be assumed constant and equal to i'i63. 
In Table I. are given the values of £ and the moduli of C, D, *E, 8E 

and P for various values of hY. The values of SE = E-—£P 3 

1163 

are given in order to show the accuracy of the approximation 
made. 


Table I. 


hV . 

p. 

f. 

c. 

D. 

E. 

5 E. 

0*0 

*(3000 

1‘585 

•0000 

•OOOO 

•OOOO 

'OOO 

o-i 

•0566 

1 573 

•0884 

•0050 

*0002 

•OOO 

0*2 

•1124 

1*561 

*1755 

-0197 

*0017 

•000 

0*3 

•1680 

1'548 

•2600 

*0437 

•0057 

- *OOI 

0*4 

•2227 

r ‘53 1 

*3409 

*0759 

•OI34 

•001 

o*s 

•2762 

1*510 

•4171 

•1152 

•0253 

•002 

0*6 

•3284 

1*486 

•4879 

•1602 

•0420 

•003 

07 

•3789 

1’459 

*5528 

•2094 


•005 

o # 8 

•427s 

* "43* 

•6115 

*2614 

•O9OI 

•006 

0*9 

•4739 

1*401 

•6640 

•3147 


*007 

i # o 


1*371 

•7103 

*3679 

•1558 

•008 

ri 

•5595 

1*342 

•7507 

-4200 


•008 

1 *2 

•5984 

I ’3 I 3 

•7856 

•4701 

•2337 

•008 

I# 3 

•6345 

1*285 

•8156 

•5175 


*007 

I# 4 

•6678 

1*259 

•8410 

*5617 

•3172 

•005 


'6984 

i'235 

•8626 

•6024 

*3590 

- ‘003 

1*6 

•7262 

1*213 

-8807 

•6396 

•4OOO 

•O'JO 

17 

7514 

1*192 

•8959 

*6732 


+ *004 

i*8 

7742 

1*174 

•9087 

7035 

*4774 

•009 

1*9 

•7946 

i*i57 

•9194 

•7305 


*014 

2*0 

•8128 

1-142 

•9284 

7546 

•5468 

+ ‘019 


With this substitution, (io) becomes 


2L = (1 + a) £ x P x + (1 - a) £ 2 P 2 ; etc. 

Now let (1 + a)£ 1 = (i + fi)Jc 1 

(!-«)&=(■ -mi 

so that + a)£ 1 + |(i - a)£ 2 . 

We obtain finally 

2 L/* = (i+/3)P 1 + (i-/3)P 2 . 
2M/A = (i+j8)P 1 a + (i-j8)P 2 a 

2 N 7 /* = (i+/3)P 1 3 + (i-^)P 2 3 


. (12) 

• ( r 3 > 


• (14) 

• (15) 

. (16) 
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It will be seen from (13) that A; is a weighted mean between 

and £ 2 , the weights being proportional to the number of stars in 
the two drifts. Now looking at Table I. we see that £ z and £ 2 vary 
within fairly narrow limits; actually the values of hv 1 and hv 2 in 
the different parts of the sky are such that, except in a few unusual 
cases, k will lie between 1*35 and 1*45. For a first approximation, 
we may either calculate a provisional value of k from the roughly 
known speeds of the drifts, or we may assume for it the value 
i‘40, which will certainly be not far from the truth. A second 
approximation may afterwards be made if desired, using the value 
of k deduced from the results of the first approximation; in 
practice, however, this is usually quite unnecessary. 

Thus k may be assumed to be known, and the solution of (14), 
(15), and (16) is then easy. We find 

IW-tMHi-rtP.-Pffi + P,) . (17) 

therefore N-)*-L M p R 

M&-L 2 12 

whence solving (14) and (19) for P x and P 2 

■P.-K + d-K)/*! 

P2 “ k -S - k )/^ ' 

and substituting these values in (15) we easily find 

L - M 

' V(L-K4) 2 + (M*-L ! ) 

We also find from (14) and (18) 

* 

p --r c ( L+ v / 7r|' /MFrE:2 )i 

Equations (21) and (22) constitute the solution of the problem. 

The argument of P z is 6 V and hv l can be found from the 
modulus by means of Table I. £ z and £ 2 can be taken from the 
same table, and then a is found from the formula 

' I - a 1 — )8 £ z 

I + a i 4 - fi £ 2 


• (i9) 


. (20I 


(21) 


(22) 
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In the practical application of these formulae a difficulty arises 
from the fact that, whereas f 3 is by definition a real quantity, the 
right-hand side of (21), which is derived from observation, is in 
general a complex quantity. 

The reason is this: equations (14), (15), and (16) practically 
constitute six equations, for each of them has a real and an 
imaginary part, which must be satisfied separately. But there are 
only five unknowns. In the above solution the equations have 
been reconciled by introducing a sixth fictitious unknown, viz. the 
imaginary part of a. If, however, a and (3 are to be purely real, 
the six equations cannot be exactly satisfied. Now it is clear that, 
just as we were, justified in rejecting equations derived from equat¬ 
ing the higher harmonics, we should attach more weight to equations 
(14) and (15) than to (16). It seems reasonable to satisfy the 
former exactly, and to throw the discordance entirely on (16). In 
other words, the observed and theoretical distributions will be mad6 
to agree exactly as regards coefficients of sin 0, cos $, sin 26, cos 26, 
in the corresponding Fourier series, the coefficients of sin 3 6 and 
cos 3$ will be made to agree as well as possible, and all higher 
harmonics will, as before, be entirely disregarded. If this rule 
seems rather arbitrary, it must be remembered that the whole 
question is one of weighting the equations of condition, that the 
ideal system of weighting is impracticable, and that we are seeking 
a practicable system which shall not unduly waste the material 
afforded by observation. 

As (16) is not to be exactly satisfied, our solution must be chosen 
to make the sum of the squares of the two residuals (of the real 
and imaginary parts of the equation) a minimum. This is equiva¬ 
lent to making a minimum the modulus of the residual of the 
complex equation. Or we may satisfy the equation exactly by 
replacing N by N 4 - 8N, where | SN | is to be a minimum. 

Following this change through the solution previously given, 
let the alteration in N change L - KA to (L - K&) 0 ; we, see from 
(19) that if | SN | is a minimum, 


| (L - K&) - (L - K &) 0 1 will be a minimum. 
But (21) may be written 



therefore if (3 is real, we must choose (L — KA) 0 so that its argument 
is the same as, or differs by 180° from, that of - L 2 . 

(The alternative arg (L - K&) Q = arg jWc - L 2 ± 90° leads to a 
value of (3 greater than unity, which must evidently be excluded.) 

The conditions are satisfied if (L - K/t) 0 is the projection of the 
vector L — JK k (in an Argand diagram) on the direction of the 
vector ± JM-k - L 2 ; for this makes | L - K& - (L - K&) 0 | a 
minimum, subject to (3 being real. 
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Therefore 


| (L - K &) 0 1 = | L - Kk | cos {arg (L - K&) - arg jMk - L 2 } 


and 



1 M& - L 2 

l(L-K*) 0 


2 


( 23 ) 

( 2 4 ) 


and Pj and P 2 are given by (22) without modification. 

The sign of (3 is positive or negative according as (L - K&) 0 is 
of the same or opposite sign to Juic - L 2 . Either square root of 
M& - L 2 may be taken, but the same root must be adhered to 
throughout. It is, however, desirable to choose the root whose 
argument agrees roughly with the direction of motion of Drift I 
relative to Drift II, otherwise the suffixes of the constants will not 
agree with the usual designations of the drifts. 


Example. 

The following example shows the practical application of the 
method. It refers to Region C (R.A. 2 h -6 h , N.P.D. 20°-52°) of 
the Groombridge proper motions. 

From the observed proper motions I find 


Total number of stars n — 512. 

2 cos 0=+73 - o Sain #=+2i9'o 

2 cos 20= — 42’3 2 sin 26— + 124/1 

5 cos 2)0= - 98’45 S sin 3^=+ 23*1 

hence, dividing by n, 

L — + "T43 "b ’428 L 
M= — ‘083 + *242 6 
N — — *192 + ‘045 t 

and -yH = — '223 + ’052 1 (y=i'i63•). 

Assume provisionally k = 1 ’40, we find 


LM = - ‘1967 + '07371= *2100 expi* i59°*5 
M&-L 2 = 4 - '0466 + '21641= '2214 expi 77°'9 
therefore by division (equation (19)) 

2K= ‘948 expi 8i°'6 

K — + *069 + '469 6 

therefore L - = + '046 -‘2291 

= '234 expi 28 i° - 4 ; 

also Jm. k - L 2 = '470 expi 39°‘o. 


If we continue the solution witli L - K k unmodified, equation 
(21) gives 


( 3 = - '300 - '442 l. 


* Expi x = e lX . 
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and 


But modifying it in accordance with (23) 

(L - K&) 0 = '234 cos (28 i°' 4 - 39°*o) expi 39°‘o 

= - '1083 expi 39°’o 

(L-.K*) 0 - _ 
whence [3= - 


4 - 53 . 


•A+( 4 ' 53 ) s ' 2 ' 5 ' 

This value of (S must now be substituted in (22), 


v 7 ^- 1 ' 244 - 


J 


1 -P 


JWc - L 2 = ’610 expi 39°*o 


1 +/3 

-- *474 + '384 1 
and L = *143 + *4281. 

Sum = M > 1 = *6i7 4-*8i2 t, 

= 1*020 expi 52**8. 

Therefore Pi = * 7 2 9 ex pi 5 2 ° *8- 

Similarly we find P 2 = *i74expi 
The directions of the two drifts are accordingly 

0 x = 52°*8 0 2 =i32°*7. 

To find their velocities we make use of Table I., which gives 
the corresponding values of [ P | and hv, 


hv 1 — i*6i 


also from Table I., 


but 


£1=1-211 

1 -a- = 1 ~/3 Ci 
1 + a i 4- (i £ 2 


hv 2 = 0*31; 

4 = 1 ’547 


= 1*548 x 


x *211 
1 '547 


= 1*212, 


whence a— — *096 

^(i+a) = *452 J(x-a) = *548. 

Thus 45*2 per cent, of the stars belong to Drift I and 54*8 per 
cent, to Drift II. 

Finally, we can check the provisional value of Jc, 

& = i(i+a)£i + !( r - a K 2 =I '395> 
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agreeing with the adopted value 1 '40, so that no further approxi¬ 
mation is needed. 

The results obtained by this method agree closely with those 
previously obtained by the method of trial and error (“ Systematic 
Motions,” p. 50). I give below a comparison of the constants deter¬ 
mined in the two ways for four of the G-roombridge Regions; it 
affords some sort of indication of the reliability of both methods. 




Analytical Method. 

K 1 

Trial and Error. 

A 



Drift I. 

Drift II. 

Drift I. 

Drift IL 


/ hv 

1-38 

0-31 

I’65 

0-30 

Region A 


o° 

15 7° 

o° 

I 55° 


in * 

52-5 

47 ‘S 

487 

5i'5 


( hv 

* *47 

0*48 

i'55 

o'45 

Region B 


12° 

13P 

IO° 

130° 


( n 

49*5 

5 ° ‘5 

49'5 

5°-5 


f hv 

r6i 

0-31 

1'65 

0*30 

Region C 

. ] e 

53° 

r 33 0 

55° 

125 0 


( n 

45 - 2 

54-8 

4**5 

58-5 


f hv 

o-86 

077 

1 *20 

o’45 

Region F 


226° 

75° 

225° 

8o° 


( n 

63-0 

37‘0 

465 

53*5 


* Per cent. 


The weakest point in the solution is exemplified in the case of 
Region F; the determination of a (or ft) is not very satisfactory. 
The difference in the two solutions in the case of Region F depends 
entirely on the division of the stars between the two drifts; if we 
had adopted the same value of a in the two cases, the remaining 
four constants found from (22) would have agreed almost exactly 
with those found by trial and error. The weakness of the deter¬ 
mination of ft can be shown analytically. 

Differentiating (21), * 

-^M(M*-L3) + (L-KA)2}-f 


dft 
dj sr 


Now | P x - P 2 | is generally about 07, but may be less, hence 


is generally about 1 o, but may be more. 


Thus in a region containing 500 stars, a change of 2*5 in 
2 cos 3 6 or 2 sin 36 would produce a change in N of ^005, and 
the corresponding change in ft or a might be ± -05. 

The question arises whether the weakness of the determination 
of ft is a defect of the method of analysis or is necessarily involved 
in the nature of the observational data. I think there, is little 
doubt that the latter alternative is correct; I have examined 
various other methods of determining ft, but all are rather insensi¬ 
tive. In the case of Region F, both the solutions given above are 
found to agree with the observed distributions almost equally /welL 
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Thus it seems likely that even in an ideal solution we should have 
one equation very much weaker than the other four (the weakness 
may be more pronounced for some regions of the sky than for 
others). As the five constants depend on one another and on this 
fifth equation, they may all share in the uncertainty. 

One or two considerations help to avoid this difficulty to some 
extent. We may be content to assume a=o, i.e. that the stars are 
evenly divided between the two drifts; all evidence seems to in¬ 
dicate that this is approximately true, and it is conceivable that 
there may be some physical reason for it. Or, instead of adopting 
the value of a found for thh particular region, we may adopt a 
mean found from all the regions discussed; this will have a much 
smaller probable error. But the most fortunate circumstance is 
that we mag determine the relative motion of the two drifts almost 
independently of a. 

Equation (18) gives 

Pi - P 2 = -£-( 1 - P 2 )-* jMk - U . 

tc 

If p lies anywhere between -0*3 and + o - 3, (1 may be 

put equal to *98 with an error certainly less than 3 per cent.; but 
this range of p includes all values likely to occur. Only in the 
ease of a very great disparity in the distribution between the two 
drifts could a value outside these limits occur. Thus Pj P 2 is 
nearly independent of p or a. P is a sufficiently nearly linear 
function of hY for h(Y 1 - V 2 ), to be also nearly independent of a. 
Thus although adopting a = o may lead to some error in the 
determinations of hY l and hY 2 , the error will nearly be eliminated 
.from the determination h(Y 1 - V 2 ). This relative motion of the 
two drifts is the quantity which most interests us, especially as a 
systematic error in the proper motions does not affect its determina¬ 
tion so adversely as it affects the determinations of hY 1 and hY 2 . 

Systematic Motions of Zodiacal Stars. 

I have applied the theory given above to the discussion of the 
proper motions of the zodiacal stars. The proper motions were 
taken from the Catalogue of Zodiacal Stars, Astronomical Papers 
of the American Ephemeris, vol. viii., part iii. Excluding the 
stars of the Pleiades, this contains 1533 proper motions. I divided 
the zodiac into sixteen regions, each extending 22^° in longitude 
by about 16° in latitude; these are denoted successively by la, 

Ila,. Villa, lb , .VIII6, the centre of region la being 

at the first point of Aries. As regions la and 16 are diametrically 
opposite to one another, the observed motions are in parallel planes, 
and the two regions may be treated together; similarly, the other 
regions can be treated together in pairs. Thus the number of 
regions is virtually reduced to eight, each containing from 150 to 
250 stars. 

The distribution of the proper motions as regards direction in 
the eight regions is shown in Table II. Opposite 6 = 0° in the 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 






_3_8 8 S *_*_8 9 1 *_S¥HNW8_0 61 


June 1908. Mathematical Theory of Two Star-drifts, etc. 599 

first column is given the number of stars whose observed motions 
are in directions between # = 355° end # = 5°, and so on. The 
numbers have not been smoothed. For stars in la, Ila, etc., 
0 =o° is in the direction of increasing R.A., and 6 = 90° in the 
direction of increasing Dec. For stars in 16 , 116 , etc., the reckon¬ 
ing of 0 agrees with that in the opposite parallel planes la, Ila, etc. 

Table II. 

Distribution of the Proper Motions In Direction. 


Regions. 


1 0 

I. 

II. 

III. 

IV. 

V. 

VI. 

VII. 

VIII? 

'■ o ° 

16 

14 

9 

5 

7 

2 

1 

2 

10 

11 

8 

3 

8 

0 

O 

5 

4 

20 

3 

■ 3 

3 

2 

1 

I 

2 

2 

30 

2 

6 

O 

2 

3 

4 

2 

I 

40 

I 

2 

2 

0 

1 

2 

0 

2 

50 

4 

3 

O 

2 

3 

2 

3 

O 

60 

2 

I 

O 

1 

0 

1 

1 

O 

70 

3 

2 

O 

0 

0 

0 

2 

2 

80 

2 

I 

O 

1 

1 

0 

0 

2 

90 

3 

3 

I 

1 

1 

0 

1 

I 

IOO 

1 

1 

I 

0 

1 

1 

0 

O 

no 

1 

1 

O 

2 

’ 1 

1 

2 

O 

120 

0 

2 

I 

1 

0 

4 

2 

3 

130 

0 

1 

I 

2 

3 

2 

2 

4 

140 

1 

0 

I 

2 

0 

3 

5 

4 

150 

2 

1 

O 

1 

4 

1 

3 

5 

l60 

3 

2 

I 

2 

5 ' 

6 

4 

3 

170 

1 

2 

3 

2 

3 

2 

9 

11 

l8o 

0 

X 

2 

2 

6 

10 

8 

13 

I9O 

6 

4 

1 

1 

6 

7 

18 

M 

200 

6 

8 

4 

3 

8 

* 9 

12 

12 

210 

5 

5 

6 

3 

5 

9 

12 

*4 

220 

6 

4 

3 

5 

11 

j6 

14 

8 

23O 

10 

5 

6 

8 

11 

20 

10 

$ 

24O 

2 

7 

2 

4 

13 

14 

16 

4 

250 

2 

6 

3 

10 

9 

17 

13 

4 

260 

1 

,3 

6 

7 

9 

9 

9 

5 

27O 

6 

5 

2 

7 

21 • 

12 

3 

8 

280 

3 

6 

9 

12 

17 

10 

6 

4 

29O 

4 

7 

13 

20 

10 

8 

10 

4 

3 °° 

6 

7 

11 

18 

16 

8 

0 

10 

310 

12 

7 

8 

13 

13 

6 

7 

3 

320 

17 

11 

'15 

27 

13 

6 

2 

2 

33 ° 

14 

4 

12 

18 

8 

4 

2 

5 

340 

9 

20 

18 

23 

6 

1 

0 

3 

35 ° 

10 

17 

9 

25, 

5 

0 

8 

4 

Total 

175 

180 

156 

240 

221 

198 

194 

169 
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.Mr. A. S. Eddington , On the 


LX VIII. 8, 


It will be seen, by looking down the columns of the Table 
that the two streams are plainly marked in Eegions I, II, III, and 
VIII. In the other four regions their directions (projected on the 
sky) are inclined at an acute angle, and the existence of the two 
streams is rather concealed. The success of the analysis in these 
four regions is on that account especially interesting. It may be 
noticed that this belt of the sky is not so favourable as the Groom- 
bridge region for showing conspicuously the separation into two 
streams ; the centre of the latter region lies between the two apices, 
so that in it the streams are in nearly opposite directions. 

Although, in the main, pairs of regions such as I a and I& were 
treated together, I thought it safer to calculate JMk - L 2 and L 
for I a and I& separately, and to take the mean afterwards. This 
was in order to avoid the possible effects of systematic error, by 
ensuring that the difference of motion of the two drifts found from 
the observations was a difference of motion of intermingled systems 
of stars, and not the difference in the apparent motion of stars in 
la from that of stars in I&. Actually, however, the precaution 
might have been omitted. I found that in every region very nearly 
the same result was obtained whether the two halves were treated 
separately or together. 

For the solutions I used entirely the equations— 




I did not calculate N, or attempt to find yS from the observations. 
The final results given below depend, therefore, on the assumption 
(based on previous experience) that the stars are evenly distributed 
between the two drifts; but calculations are given which show to 
what extent the results obtained would need to be modified if the 
assumption is incorrect. 

As a preliminary I made two solutions, (a) assuming 
/ 1 ~ o = 1X1 an d ( b) assuming /- —|=ri; in both cases 


k was assumed to be 1*40. These correspond to assuming that 
the stars belonging to Drift I are about (a) 55 per cent., (6) 50 
per cent, of the whole. 

The combined results from all the regions were— 


Velocity of Drift I { ^ 1 

* 1 ( 6 ) 174 


1 '61 towards Latitude 




Velocity of Drift II -f [ a } 0 ^ 
J 1(6) o-6o 


302° 

304 ° 

I0 9°'4 

io9°*6 


Velocity of Drift I ((a) 1*84 
relative to Drift II 1 ( 6 ) i’9o 


relative to Drift II 1 ( 6 ) i’9o 


of Drift I / (a) 1 '84 


The unit of velocity is i/h. 
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The mutual relative velocity of the two drifts is thus determined 
nearly independently of the assumed division of the stars between 
them, a result which has already been arrived at theoretically.' 

For a final solution, using ( b) as a first approximation, I cal¬ 
culated for each region values of k and / 3 , assuming that the stars 
were equally divided between the drifts. 

The values of k for the Regions I, II, ... . VIII were 
respectively— 

. i*35. 1 'S 5 , i *3<5, i *4°i i* 44 , 1 ‘44, i‘43> I*3 8 , 



1 "x2, 1 ’12, i*ii, 1*09, 1*07, 1*07, 1*08, i'io. 

Performing the analysis with these values, the constants of the 
drifts in the eight regions were found as follows:— 


Region. 

Longitude of Centre. 


Drift I. 

Drift II. 

_X_ 

I 

0° 

, 180 0 

r 

A«i. 

i *57 

eT" 

345°‘4 

hv 2 . 
0*46 

■N 

#2- 

244 °’ 8 . 

II 

22% 

, 202% 

1’46 

350 -o 

0*58 

249 *4 

III 

45 

, 225 

i’86 

330 -o 

0*63 

252 ’2 

IV 

m 

, 2 47 i 

i-68 

323 ‘5 

0-51 

271 -3 

V 

90 

, 270 

i’o8 

299 7 

070 

23I "2 

VI 

112% 

. 2924 

1’04 

245 ‘4 

0*54 

233 ‘3 

VII 

*35 

> 315 

* "35 

.208 ’8 

0*40 

275 

VIII 

1574 

» 337 £ 

1-49 

189 ’2 

°‘54 

292 ’3 


If we resolve these drift-velocities along and perpendicular to 
the plane of the ecliptic, the componerits perpendicular to the 
ecliptic derived from each region should agree. 

The determinations of these components are— 

for Drift I, - *97, - 77, - 1 -35, - 1 ’19, - *94, - 1 -oo, - ’96, - 76, Mean - *99 
for Drift II, - -30, - -43, - -52, - *51, - *54, - ’48, - -37, - ’39, Mean - -44 

The accordance of these seems decidedly good when it is 
remembered that they are derived from regions containing on an 
average less than 200 stars. 

The component drift-velocities in the ecliptic are more or less 
foreshortened according to the longitude of the region; I deter¬ 
mined the mean values by a least-squares solution. 

The final results are— 


Velocity of Drift I, 

00 

1—1 

towards Latitude 

-33 0, «, 

Longitude 

x°5°’5 

Velocity of Drift II, 

0-59 

1 

>> 

> 5 

-47°-9. 

M 

3°3° 7 

Velocity of Drift I, rela-' 
tive to Drift II, 

^ 1 *94 

9 ) 

J J 

- i6°*4, 

i ? 

I 0 9°‘3 






43 



© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 







_3_8 8 S *_*_8 9 1 *_S¥HNW8_0 61 


602 Mr. A. S. Eddington, On the lxviil 8 , 

The slight difference between this result and that of solution 
(b) represents the effect of allowing for the variation of fi and k 
from region to region, instead of adopting a mean. The comparison 
between solutions (a) and (6) serves to indicate how the results 
would be changed if the stars were not evenly divided between the 
drifts. The constants most affected by such a change are the speed 
of Drift I and the direction of Drift II; the other constants are 
nearly independent of the assumption, and are therefore more reli¬ 
ably determined.* 

From an examination of residuals, I estimate that the probable 
accidental error of the determination of hY x and hV 2 is about ± '06, 
and the probable errors of the apices of Drifts I and II are re¬ 
spectively about 2° and 6° of a great circle. 

Converting latitude and longitude into R.A. and Dec., the 
positions of the antapices may be compared with previous deter¬ 
minations as follows— 



Drift I. 


Drift II. 


Kapteyn 

. R.A. 85° 

Dec. 

-n° 

R.A.. 260° 

Dec. 

- 4 8 >o 

Dyson . 

94° 


- 7 ° 

>, 2 40° 

>) 

1 

0 

Groombridge stars 

,, 90° 

j j 

-19 0 

„ 292° 

? > 

0 

00 

1-0 

1 

Zodiacal stars 

„ 103 0 

>> 

-ii°- 

» 330° 

? > 

-6 4 ° 


The great R.A. of the antapices of both drifts found in the present 
discussion is rather hard to account for. 

The velocities of the two drifts hv x — 178, hv 2 = 0*59 are in 
excellent agreement with those found from the Groombridge stars 
7^ = 17, hv s = 0-5. 

For the velocity of one drift relative to the other, the Zodiacal 
stars give the value i'94. • From the Groombridge stars (by a least- 
squares solution from the results of the separate regions) I have 
found the value i’c)o. The determinations of the point towards 
which this relative velocity is directed (called by Professor Kapteyn 
the true vertex) are 

Kapteyn .... R.A. 91 0 Dec. + 13° 

From Groombridge stars . ,, 95° ,, + 3” 

„ Zodiacal „ . „ 109° „ + 6° 

Professor Schwarzschild’s determination of the line of symmetry 
of motion may be added. This, although based on a rather different 
theory, is directly comparable with the above. He found, 

From Groombridge stars . .R.A. 93 0 Dec. + 6° 

* It should be understood that it is not impossible to determine o from 
the observations, bub simply that when, as in the present case, a few stars 
are discussed, the value of a is liable to a greater uncertainty than some of 
the other results. Undoubtedly from the whole 1533 stars a fairly good 
mean value of a could be determined; but the task of computing it for sixteen 
regions separately, arid taking the mean, would be laborious. 
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Further, from Professor Dyson’s investigation, an R. A. of about 
92 0 for this point may be inferred. 

Thus the Right Ascension given by the Zodiacal proper motions 
is discordant as compared with the other determinations. The 
number of stars here considered is fewer, and the proper motions 
are perhaps not so well determined, but I do not think the dis¬ 
cordance can be altogether attributed to this. Nor can it be traced 
to a local anomaly, for it seems to be indicated by the proper 
motions all round the ecliptic. I have verified by calculation that 
the Regions la to Villa agree with Regions 16 to VIII6 in leading 
to this high value of the Right Ascension. 


Note on the evaluation of certain integrals required 
in the analysis. 

f 27 r 

(1) To calculate n 1 C 1 = p cos (0 — 6 f)d6, for a single drift. 

Jo' 

We may choose the initial line, along the direction of the drift, so 
that 6 1 = 0. 


7r 


j Jjr + rer2 J e ^dx | 


Then 
where t = hv l cos 0 . 

Now 1 1 + rer 2 J e~ x ‘‘dx j cos 6 is an odd function of cos 6 and 
vanishes when integrated from o to 27r. 


/’ 


Hence I p cos 6 d6 = f 

7 T 2 J 


re* 2 cos OdO 


o 

27 r 


cos a^gW^aoos^. 


Now 


f 


e"C0S2 6dQ 


2 sj 7T 

rr r> 

e 2KC0a2d d$ = e K 
J <_ 

= 27 re K J 0 (tK). 

Whence differentiating with respect to k, 

f 27 r 

2 I cos 2 $e 2K 008 20 dO = 27re K (J 0 (ck) - iJJik)). 

Jo * 

Now write k = and substitute above 

e,= jafhv ie -*** {(*PV)} 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 






_3_8 8 S *_*_8 9 1 *_S¥HNW8_0 61 


604 Mathematical Theory of Two Star-drifts, etc. lxviii. 8 , 

Tables of the Bessel Functions of an imaginary argument are 
given in Brit. Assoc. Reports, 1893 and 1896, and have been used 
in calculating Table I. 

C 27 r 

(2) To calculate WjDj = I p cos 20d0. 

J O 

We have p — f e ~^ 2 (^ 2 +vi 2 - 2 ^i rcos 0)r<#r 

7 T J O 

and identically 


= f pdO = -( d6 f e-i¥+«i 2 -2¥co e )rdr 
o 7 T J o Jo 

= Pe-^rdr f 2 V 008 w. 

7T J o •/ o 

where 2 = zhfiv^r. 

Hence T e -^+*\dr.2TrJ 0 (iz) . 

'TT J o 

But from (1) 


(1) 


(2) 


/: 


p cos 2 0dd 


nf 2 


7 T 


i 


e-w^+^rdr 


,df f” 

dz 2 J o 


e z COB 9d0 


= V ^L r e-’^+^rdr—izTrJfiz)) 

7T y o 

= ^ e -ft2(r '‘' i+1 ' l2) rc?r^27rJ 0 (iz) — i^(27rJ 0 («,2))^ 


by Bessel’s differential equation 


7 T J o 


9 —^(r^+vi 2 ) 




4 / t 4 ?;-, 2 dr 

Integrate by parts: 

+ *_ rdr(- 2h 2 r)e- h ^ + ^2- n -JMz) 

it 4 h'hfj o 


= M 4. jLr^i’_ n i_ 

1 2 h 2 vf 2 h 2 vf 


71 

(the integrated part vanishes for r = 00 and reduces to 
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at the lower limit; the part remaining to be integrated is simplified 
by means of the identity (2)), 


hence 


i‘ 


p cos 2 6 d6 ■■ 


p cos 2 6 d 0 


f 27 T 

I pd 6 

J O 


= 1 - 


1 - e 


-l&V? 


h*vf 


f 27 T 

(3) The integral I p cos 3 6d6, and corresponding integrals for 
J O 

any odd multiples of 6 , can be found by a simple extension of 

/* 27 T A 

the method employed for I p cos 6 d6 . The integrals for evea 

J o 

multiples of 6 are more troublesome to evaluate, but the method 
employed for / p cos 2dd0 always succeeds. 

J O 


Tables of the two hypergeometrical functions, F (1/6, 5/6, 2, sin 2 — y 

/ 1 \ ' 2 / 

and F( -1/6, 7/6, 2, sin 2 — j, between the limits iota equals 
go and 180 degrees. By C. J. Merfield. 

In the method of Mr. It. T. A. Innes for the determination of 
the secular perturbations,* there are two hypergeometrical functions 
to be deduced. In an appendix f to this valuable paper, tables of 
the logarithms of these functions are given with the argument 
iota for each degree | for the first quadrant. 

Tables of these functions facilitate the applicatiou of this 
method in no small degree, and it seemed desirable to extend 
them, as in many future investigations it will be found that the 
angle iota will much exceed a quadrant. Taking an example, 
Eros §—Earth, it will be noted that the modular angle theta, the 
argument to the tables of elliptical integrals, reaches the value 60“, 
corresponding to iota 138° 18', and there will be many other cases 
in which it exceeds this value. 

In the preparation of the tables here given the formulae || (Z.c., 

* “Computation of Secular Perturbations,”by R* T. A. Innes, Monthly 
Notices, vol. lxvii. 427. 

t Tables for the application of Mr. Innes’s Method, by Frank Robbins. 
l.c», 444 * ' 

t The values of these functions for iota equals o° have been omitted 
in the tabulation by Mr. Robbins. 

§ “Secular Perturbations of Eros,” by C. J. Merfield, Astr.Nachr., 4178-79, 
Band 175. 

|| The values of these functions may be deduced from series. I have given 
the coefficients of twenty terms, Astr. Nachr., 4215, Band 176, p. 246. 
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